First orbital solution for the non-thermal emitter Cyg OB2 #9 
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ABSTRACT 

After the first detection of its binary nature, the spectroscopic monitoring of the non-thermal 
radio emitter CygOB2#9 (P=2.4yrs) has continued, doubling the number of available spectra 
of the star. Since the discovery paper of 2008, a second periastron passage has occurred in 
February 2009. Using a variety of techniques, the radial velocities could be estimated and a first, 
preliminary orbital solution was derived from the He i A 5876A line. The mass ratio appears close 
to unity and the eccentricity is large, 0.7-0.75. X-ray data from 2004 and 2007 are also analyzed 
in quest of peculiarities linked to binarity. The observations reveal no large overluminosity nor 
strong hardness, but it must be noted that the high-energy data were taken after the periastron 
passage, at a time where colliding wind emission may be low. Some unusual X-ray variability is 
however detected, with a 10% fiux decrease between 2004 and 2007. To clarify their origin and 
find a more obvious signature of the wind- wind collision, additional data, taken at periastron and 
close to it, are needed. 

Subject headings: binaries: spectroscopic - stars: early-type - stars: emission-line - stars: individual 
(CygOB2#9) - X-rays: stars 



1. Introduction 

As the main responsibles for mechanical in- 
put, chemical enrichment and ionizing radiation in 
galaxies, massive stars (i.e., with masses >1OM0) 
are important objects of the stellar population. 
However, these objects are rare: in the Galaxy, 
about 300 Wolf-Rayet stars (i.e., fewer than ex- 



oplanets) are known (jvan der Huchtl l2006l) . as 
well as about 400 0-type stars (plus 700 can- 
didates that h ave so metimes be classified as O, 
2008h . 
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A natural consequence is 
that many aspects of these stars remain poorly 
known, even in the case of their major property, 
the stellar wind, for which "basic" quantities such 
as the mass- loss rates are still heavily deba ted at 
the present time (e.g. Sundqvist et al. 2010l ). 

In this context, non-thermal radio emitters as- 
sociated with massive stars form an even more 
limited group, with fewer than 40 cases known in 



iimitea grou p, witn lewer tnan 4U cases xno wn m 
our Galaxy jPe Beckeij[2007l : lBenagliall2009l ). but 
these objects can provide unique insights into the 
physics of stellar winds since both phenomena are 
intimately linked. 

Observing such non-thermal radio emission im- 
plies two pre-requisites: the presence of both a 
magnetic field and a population of relativistic elec- 
trons. Direct detection of the former is notoriously 
difficult in massive stars, which display only few 
broad lines and therefore have a weak Zeeman sig- 
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nature. However, as some of the descendants of 
massive stars (i.e., neutron stars) are clearly mag- 
netic, there was little doubt that magnetic fields 
are not totally absent in these objects. Indeed, the 
last decade saw the collection of the first evidences 
of magnetism in massive stars tha nks to spec 



tropo larimetric monitorings (e.g., iDonati et al 
2002. with even a tentative detection in a non- 



thermal emitter in lHubrig et al.1120081 ). 

The second requirement, the population of rel- 
ativistic electrons, is ultimately linked to the pres- 
ence of shocks in the ionized winds (for a re- 



De BeckedlioOTh . Acceleration then takes 



view see 

place through the first-order Fermi mechanism, 
also called "diffusive shock acceleration", where 
electrons iteratively gain energy by crossing sev- 
eral times the shock zone. Shocks are not ubiqui- 
tous in stellar winds: they can either be intrinsic 
to the winds themselves (as the line-driven mecha- 
nism producing the wind is intrinsically unstable) 
or arise in massive binaries from the collision of the 
two stellar winds. Recent theoretical modelling 
showed that only the latter can acco unt for the ob- 
served non-thermal radio emission (van Loo et al 



20051 ). Indeed, the large majority of the known 
cases of non-thermal emitters arc confirmed or sus- 
pected binaries - only 3 of the 17 WRs a nd 3 of 
the 16 O stars listed by IPe Becked (|2007t ) totally 
lack evidence for binarity. However, this can of- 
ten be explained by an inadequate monitoring, or 
sometimes simply the absence of any monitoring. 
Once an adequate observing campaign is orga- 
nized, proofs of binarity are often found, as exem- 
plified by our re cent successes in this d omain (e.g., 
CYgOB2#8A. iDe Becker et al.l I2004L and 9Sgr 
Rauw et ahlbOOsl 2010. in prep) which greatly im- 
proved the multiplicity census of non-thermal ra- 
dio emitters associated with 0-type stars. Among 
the remaining objects without evidence of multi- 
plicity, the early-type star CygOB2#9 (051) is 
clearly a target of choice. 

Located in the rich association Cyg 0B2 

(|Kn6dlseded l2000l ). CygOB2#9 was one of the 
first 0-stars shown to be a non-thermal radio- 
1981 



emitter (jAbbott et al 
this emission 



The 



presence of 
problematic for years. 



remained 

The first direct evidence of the binary nature of 
Cyg 0B2 #9 was only reported in 2008 thanks 
to a d edicated, long-ter m spectroscopic moni- 
toring ( Naze et aD 120081 ). In the same year. 



van Loo et al. I (l2008l) revealed a long-term modu- 
lation of the radio emission and derived for it a 
period of 2.355 yrs, interpreted as being associated 
with the binary orbit. This radio flux appeared 
minimum when the stellar lines were unblended 
in the optical spectrum. Up to now, however, a 
full orbital solution and a detailed modelling of 
the system are still missing. This paper aims at 
beginning to fill this gap, by tackling the first 
problem: the derivation of the orbital parameters. 

It must be recalled that Cyg 0B2 #9 is not 
an easy star to analyze. On the one hand, 
it is strongly extinguishe d {E[B — V] ^ 2.11, 
Massev fc ThompsonI 19911 ) with two consequences. 



First, the spectrum of CygOB2 7^9 displays very 
strong interstellar lines, often affecting the stel- 
lar lines, even those usually quite uncontaminated 
(e.g.. Civ A 5812). Second, the star appears faint, 
especially at blue wavelengths where most lines 
used for spectral classification are found. On the 
other hand, the stellar lines remain totally blended 
for about 80% of the orbit ; a clear doubling of 
the lines can only be seen during a few months 
of the 2.355 yrs period. This is due to the large 
eccentricity of the system but also to the width of 
the lines (FI^i/MHciA5876 3A). 

Despite these difficulties, we have continued our 
monitoring of CygOB2#9, with the hope of im- 
proving our knowledge of this system. This paper 
reports on the new data collected since January 
2008, including at orbital phases close to the pe- 
riastron passage of 2009 (the pcriastron itself was 
unobservable as Cyg 0B2 ^9 was in conjunction 
with the Sun at that time). A preliminary orbital 
solution is presented here for the first time. This 
is a necessary and crucial step towards the full 
modelling of the system and the derivation of the 
winds' parameters. 

Complementary high-energy data are also pre- 
sented here, as they directly relate to the question 
of colliding winds. At first, one could expect for 
these peculiar objects a non-thermal X-ray emis- 
sion, a direct companion to that observed in the 
radio range. However, direct and undisputed ev- 
idence of such emission is still lacking: its detec- 
tion awaits the adv ent of sensitive observa tories in 
the >10keV range |De Becker et al.ll2009[) . In the 
meantime, one could however try to find evidence 
for wind- wind collision in the thermal X-ray emis- 
sion which dominates the 0.3-10. keV range, the 
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preferred bandwidth of the current sensitive facil- 
ities (XMM-Newton, Chandra). Indeed, some bi- 
naries display strong wind-wind shocks which are 
able to produce very hot plasma, hence hard and 
bright X-rays. Moreover, in these cases, the X-ray 
emission is modu lated as the stars orb it each other 
(for a review see iGiidel fc Na^l2009^ . It is there- 
fore important to check whether the X-ray emis- 
sion of Cyg 0B2 #9 bears the signature of wind- 
wind collision. This would yield a further proof of 
the link between binarity and non-thermal radio 
emission, as well as additional constraints on the 
winds' parameters. 

The paper is organized as follows: Sections 2 
and 3 present the optical datasct and its anal- 
ysis, Section 4 details the X-ray properties of 
CygOB2#9, and Section 5 summarizes our re- 
sults. 

2. The optical dataset 

The observations obtained until ea rly January 



first order, R ^ 4000). Exposure times var- 
ied from 720 to 5400 s (generally in multiples of 
600-900 s) depending on weather conditions. All 
datasets were reduced using standard IRAF reduc - 
tion routines as outlined in lKiminki et al.l (|2007l ). 



2008 h ave been described in length in iNaze et al 



(|2008f ). and only the new data will be described 
here. These new observations were obtained at 
the Haute-Provence observatory (OHP, France) 
and the Wyoming Infrared Observatory (WIRO, 
USA); an archival dataset taken at San Pedro 
Martir (SPM, Mexico) in 2004 has been added to 
the new analysis. A journal of these observations 
is provided in Table [TJ 

At the OHP, six additional spectra in the yellow 
range and one in the blue range were obtained in 
2008 using the 1.52m telescope equipped with the 
Aurelie spectrograph (grating #3, R ^ 9000). In 
2009, the Sophie cchellc instrument {R = 35000, 
39 orders over the domain 3900 — 6900A) in- 
stalled on the 1.93m telescope observed the sys- 
tem 9 times in the high-efhciency mode. For each 
dataset, the typical exposure time was 1800-7200s 
; the observations were sometimes split into sev- 
eral individual exposures which were finally com- 
bined if taken within 1-15 days. The data were 
first reduced in a standard way, smoothed by a 
moving box average and finally normalized. Note 
that data from August 10 and 13 2009 were unus- 
able due to a technical problem. 

In Wyoming, the observations were obtained 
with the 2.3m telescope equipped with the WIRO- 
Longslit spectrograph (1800 1 mm~^ grating in 



One archival cchelle spectrum, taken with the 
2.1m telescope of SPM equipped with the Espresso 
spectrograph {R = 18000, 27 orders in over the 
range 3780-6950A), was also made available and 
added to our dataset. These data were reduced in 
a standard way using MIDAS. 

To improve the wavelength calibration, we took 
advantage of the high reddening and used several 
narrow, well-marked Diffuse Interstellar Bands 
(DIBs) close to major spectral lines. Their veloc- 
ity shifts relative to a chosen reference dataset (the 
echelle spectra taken by Sophie in October 2007) 
were measured using a cross-correlation method. 
The measured radial velocities (RVs) of the stellar 
lines (see below) were then corrected by the shifts 
derived from a close DIB (e.g., DIBA5780A for 
He I A 5876 A). This ensures that the wavelength 
calibration is correct to within 5-10 km s~^ in the 
worst case. 

In total, our dataset now comprises 20 cchelle 
spectra (Sophie at OHP, AFOSC at Asiago, 
Espresso at SPM), 8 long-slit spectra at red wave- 
lengths (Aurelie at OHP, Asiago), 10 yellow spec- 
tra (Aurelie at OHP), 1 blue spectrum (Aurelie at 
OHP), and 4 yellow-to-red, low-resolution spectra 
(WIRO). These data were taken between 2003 and 
2009, with a more intense monitoring since 2006. 

3. Towards a first orbital solution 
3.1. Radial velocities 

To derive an orbital solution, it is necessary to 
(1) choose adequate stellar lines and (2) secure 
good estimates of the radial velocities (RVs). 

Choosing reliable stellar lines is no easy task for 
CygOB2#9. With strong interstellar lines and 
low signal-to-noise ratios, many lines have to be 
discarded. While the splitting is seen near peri- 
astron for He i A 4471 and He ii A 4542, the strong 
noise in the blue domain prevents us from deriving 
good RV estimates from these lines throughout the 
whole orbit. Metal lines, such as C iv AA 5801,5812 
or O III A 5592, are strongly contaminated by inter- 
stellar lines, rendering their measurement difficult. 
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In addition, the two stars of the system appear so 
similar in spectral type that there is no obvious, 
sufficiently strong stellar line belonging to only 
one component of the pair. Also, the data do not 
always cover the same wavelength range and the 
chosen line should indeed belong to the most fre- 
quently recorded domain. Our best choice for RV 
determination is thus He i A 5876, which is simulta- 
neously strong and free from interstellar contami- 
nation. However, there is a caveat to this choice: 
this line is sometimes polluted by emission in some 
extreme 0-type stars. In our dataset, there is no 
emission above the continuum level, but it is dif- 
ficult to exclude with 100% confidence the pres- 
ence of faint emission which would slightly "fill" 
the photospheric absorptions. Our results are thus 
clearly preliminary. 

The stellar lines have first been fitted by Gaus- 
sians (usually one, with two only for the 2-3 clearly 
unblended spectra) . The last column of Table [5] 
gives the result of a single Gaussian fit to the 
blended data (i.e., 0.1 < </> < 0.9). However, 
as the li ne splitting has clearly been detected 
(jNaze et al. 2008) , at least during some part of the 
orbit, we tried to improve this RV determination 
by using the correlation method (TODCOR, sec 
Zucker fc Mazch.l994i), the disentangling m ethod 
(|Gonzalez fc Levatoll2006l : iMahv et al.ll2010f ). and 
the x^-mapping method. Used on a sole line 
(He I A 5876), the first two methods give poor re- 
sults. As TODCOR involves the convolution be- 
tween a typical line shape (i.e., a Gaussian in our 
case) and the spectrum (reduced here to one single 
line) , the peak of the cross-correlation function is 
broad, increasing the error bars on the RVs to un- 
acceptable levels. Even at maximum separation, 
TODCOR does not always yield results within 
lOkms"^ of the simple two-Gaussian fitting, and 
it was thus discarded. Disentangling the Hei 
line also proved unreliable, the main reason be- 
ing that due to the orbital configuration, the lines 
are only partially separated when the primary is 
blueshifted and the secondary redshifted. The re- 
verse situation, redshifted primary and blueshifted 
secondary, is never observed (Fig. [1] left). There- 
fore, the only solution for fitting two components 
in He I A 5876 was first to fit each component at 
maximum separations using Gaussians and then 
to shift these two Gaussians, keeping their shape 
(width and depth) constant, to find the minimum 



for each spectrum (Fig. [U right). The RVs de- 
termined with this method arc listed in the third 
and fourth columns of Table [2] When available 
and close to maximum separation, results from 
this method are similar for Hen A 6683. Fig. [5] 
shows the evolution of the RVs, the two periastron 
passages can be clearly seen, as well as the relative 
constancy of the RVs during most of the period. 
When the two components are blended, the RVs 
are still rather uncertain, as the separation is then 
much smaller than the line width. 

We have detailed above the numerous caveats 
concerning our measurements and our data (blend- 
ing during most of the orbit, possible slight con- 
tamination by emission, x^-mapping vs single 
Gaussian fit). Caution should thus apply, but it 
should also be stressed that these observations and 
measurements represent the best dataset available 
for CygOB2#9. Deriving an orbital solution us- 
ing these RVs was thus attempted. 

3.2. Orbital solution 

SBl solutions were calculated in several steps. 
First, the best- fit period was select ed using the 
Fourier method of Heck et al. I (Il985i see also re- 
marks in iGosset et al.l 120011 ) . A polynomial fit 



of the folded RV curve was then calculated and 
an approached orbital solution was derived from 
it using the best-fit result amongst so l utions cal- 
culated using methods by IWolfe et al.l (Il967l) and 
Lehmann-Filhed ( 1894 ). The former method relies 
on the derivation of the Fourier expansion (lim- 
ited to 2 terms) of the observed RV curve fol- 
lowed by an identification of the coefficients with 
those issued from a series expansion of the the- 
oretical RV curve corresponding to a Kepler or- 
bit ; the latter method uses the amplitude and 
various integrations under well defined parts of 
the RV curve to get a first estimate of the or- 
bital parameters. From that starting point, a re- 
fined, final solution was found using a Levenberg- 
Marquardt minimization on a method adapted t o 
eccentric binaries (e > 0.03) bv lSchlesinger (|l910f ). 
A check of the correctness of the orbital solu- 
tion was made based on the recent algorithm of 
Zechmeister fc Kiirster ()2009l) which, by perform- 



ing the fit in terms of the true anomaly, reduces 
the non-linearity of the RV curve to 3 free pa- 
rameters instead of 6 (the other 3 free param- 
eters then being readily derived from the peri- 
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odogram coefficients) ; only few iterations of the 
Levenberg-Marquardt method are then needed to 
reach the final solution, which, in our datascts, al- 
ways agrees with the result of the first derivation. 
Errors were estimated using the diagonal of the 
variance-covariance matrix. The SBl fitting was 
performed on the RVs of the primary, on the RVs 
of the secondary, and on the velocity difference 
between primary and secondary. 

SB2 solutions relied on the main idea of the 
Liege Orbital Solution Package (lospB In this 
method, the secondary and primary velocities are 
converted into an equivalent SBl dataset, using a 
linear orthogonal regression fit between the veloc- 
ities of the two components. After this transfor- 
mation, the equivalent SBl dataset was fitted as 
described above. 

In both cas es, the period was all owed to vary 
slightly, since Ivan Loo et al.l (|2008[ ) had a rela- 
tively large uncertainty (P = 2.355 ± 0.015 yrs) 
and our observations were taken 10 to 13 cycles 
after the radio data. Indeed, the phase of the Oct. 
2006 and Feb. 2009 periastron passages seem to 
occur at (j) = 0.95 rather than at 0.0 (Fig. [2|). 
However, since we cover only two events of maxi- 
mum separations, the formal improvement on the 
period error from our sole dataset is not very large. 

Table |3] gives the derived orbital parameters 
and their associated error bars while Fig. [3] graph- 
ically shows the results of the best SB2 solutions. 
Note that smaller weights were given to the lower- 
quality data (0.3 for WIRO and 0.6 for Asiago and 
SPM, 1 otherwise). 

The best orbital solutions rely on our best es- 
timates of the RVs, i.e., those derived from the 
X^-mapping. In each case, the best-fit period is 
slightly revised downwards but agrees well, within 
the errors, with the radio determination. As it 
takes into account both components, the SBl so- 
lution calculated on the RV difference agrees best 
with the results of the SB2 solution, though the 
parameters derived for individual SBl solutions 
are never at 3-a from the results of the SB2 so- 
lution. The mass ratio is close to unity, confirm- 
ing the si milarity of the two stars of Cyg 0B2 ^9 
found by iNaze et al ( 2008 ) who proposed spec- 



tral types of 05-1-06-7. Us ing the typical masse s 
of such stars as quoted in iMartins et al ] (|2005h . 
the masses derived from the x^-mapping imply an 
inclination of about 45-50° if both components 
are supergiants, or 55-60° if both components still 
belong to the main sequence (which is unlikely in 
view of the combined 051 spectral type) . The ec- 
centricity is large, as previously suspected, with 
a value of about 0.7-0.75. The velocities of the 
center-of-mass are quite different for both stars, as 
could be expected from the fact that the secondary 
lines are never seen on the blue side. Indeed, our 
RV curves do not cross (see Fig. |3l left) but, as 
already mentioned, the RV determination when 
lines are totally blended is difficult and a slight 
crossing (by e.g., lOkms"^) of the true RV curves 
can therefore not be totally excluded. In any case, 
it would not change the fact that the center-of- 
mass velocity of the primary is blueshifted com- 
pared to that of the secondary. This is most proba- 
bly related to the fact that the wind of the primary 
is stronger than that of the secondary, leading to 
the formation of the stellar lines of the most ex- 
treme star (the primary) not at the photosphere 
but farther, in the wind itself. 

As the RVs are uncertain when the lines are 
blended, i.e., for 80% of the orbit, we checked 
our results by calculating orbital solutions using 
the RVs from ^^"ni&PPffig only when the lines 
are unblended, i.e., near periastron, and the re- 
sults of a single Gaussian fit otherwise (i.e., for 
0.1 < (p < 0.9). The weights were halved when 
the RVs of the primary and secondary are sup- 
posed identical. These solutions are given in the 
last columns of Tableland shown in Fig. [3] (right 
panel) ; they arc identified by a "prime" sign. 
As could be expected, the difference between the 
center-of-mass velocities is reduced and the semi- 
amplitudes arc slightly enlarged. It must however 
be noted that the orbital parameters are similar, 
within the 1-a error bars (except for the center- 
of-mass velocities, which agree within 2-ct). The 
largest differences are seen in the physical param- 
eters, as the larger amplitudes naturally result 
in larger masses and semi-axes, hence suggesting 
larger inclinations (^ 70° in the case of super- 
giants, ^ 90° for main-sequence objects). 



^Tho LOSP package and a preprint describing it (Sana 
& Gossot, A&A, submitted) can be downloaded from 
[http : / / staff, science . uva.nl /~hsana/losp.html[ 
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4. The X-ray emission from Cyg OB2 #9 

X-ray emission from CygOB2 was discovered 
serendipi tously when EINSTE IN was pointed at 



CygX-3 (jHarnden et al. I I1979I) . This was actually 



the first report of X-ray emission from OB stars. 
Since then, the region has been observed several 
times by X-ray observatories. The most recent 
exposures have been obtained thanks to a moni- 
toring of Cyg 0B2 # 8a with XMM-Ne wton (see 
preliminary results in lRauw et al.ll2005l) . 



4.1. The XMM-Newton dataset 

In total, XMM-Newton provided six pointings 
centered on CygOB2 #8a (Obsld=20045, 50511, 
PI G. Rauw). The first four datasets, separated by 
10 days each and with a duration of 20 ks, were ob- 
tained in October - November 2004; the last two 
datasets, of length 30 ks, were taken three years 
later, around 2007 May 1 (Table S]). All exposures 
were obtained with the same EPIC configuration 
(full frame, medium filter). The reduction pro- 
ce ss of the first four datasets is explained in detail 
in iDe Becker et al. I (I2OO6I) . The last two were re- 
duced in a similar manner and we only repeat here 
the most important information. 

The raw data were processed with SAS ver- 
sion 6.0 package. Some bad time intervals affected 
by high background events (so-called soft-proton 
flares) were rejected. A few stray-light features 
(due to singly reflected photons) from Cyg X-3 are 
visible in the lower right corner of the images. 
However, they do not affect the most interesting 
part of the field of view. The Cyg 0B2 #9 data 
were extracted within a circle (of radius 40" for 
the first four datasets, of radius 30" for MOS and 
23.25" for pn in the last two observations) cen- 
tered on the source, whereas the background was 
extracted from a nearby 50" x 20" devoid of X- 
ray sources. The EPIC spectra were analysed with 
the xspecvll.2 software. Note that the EPIC-pn 
data of Cyg 0B2 #9 in Obs. #4 and #6 arc not 
available due to the source falling partially or to- 
tally in a CCD gap. 

4.2. X-ray properties 

In the radio range, Cyg 0B2 #9 displays a clear 
signature of non-thermal emission. Indeed, con- 
sidering its quite long period, Cyg 0B2 #9 fits in 
the so-called standard scheme for colliding-wind 



massive binaries accelerating particles, and pro- 
ducing sy nchrotron radiation in the radio domain 
(jPe Beck er 2007). However, its X-ray emission, 
revealed by our XMM-Newton observations, ap- 
pears clearly thermal in nature, as several X-ray 
lines are detected in the spectra (e.g., the iron line 
at 6.7keV, see Fig. |4]). In fact, this strong thermal 
X-ray emission could easily hide, in the 1-10. keV 
range, a faint putative non-thermal X-ray compo- 
nent due to inverse Compton scattering. To be 
unveiled, such an X-ray emission should be inves- 
tigated using observatories with a large sensitivity 
above 10 keV, where the thermal emission becomes 



negligible (|De Becker et al.ll2009l ) 



To obtain a good fit to the observed XMM- 
Newton spectra, the combination of two hot, 
optically-thin plasma was needed, and the fitted 
models are thus of the type wabsint x a6s„ind x 
{mekali + mekal2)- The first absorption corre- 
sponds to the interstellar one, fixed to Nh.ism = 
1.15 X 10^^ cm"^. This value was derived from the 



reddening of lMassev fc ThompsonI (Il99ll) and the 
gas-to-dust ratio of Bohlin et al. ( 1978[ ). To allow 
for additional, circumstellar absorption, a second 
absorbing component was added and allowed to 
vary. As the circumstellar material actually is an 
ionized wind, we used the dedicated opacity ta- 



bles fr om the wind absorption mo del oflNaze et al 



( 2004 ). as for CygnusOB#8a in iDe Becker et al 



20061 ). Table [5] lists the best- fit models obtained 
for the six pointings. Note that f^^''' corresponds 
to the dereddened flux, i.e., the flux corrected for 
the interstellar absorption. 

The overall luminosity is rather typical of 
0-stars. Us i ng th e bolomctric correction of 
Martins et al. I (I2OO5I) for an 051 star (the com- 



bined type of Cyg 0B2 #9) as well as the V mag- 
nitude and reddening from iMassev fc ThompsonI 
(jl99l[ ). the log(Lx/iBOL) is found to be —6.3, 
close to the typical value of this ratio for 0-type 
stars (—6.45 w ith a dispe rsion of 0.51 dex in the 
2XMM survev lNazellioogI) . The average temper- 
ature (< kT >= '^{kTi X norrrii) / '^{normi)) 
derived from the fits is 1.2 keV, which is slightly 
high compared to the average temperatures de- 
rived in the 2XMM (where 83% of the objects have 
< kT > below IkeV). In fact, when accounting 
for the severe ISM absorption, the spectral shape 
appears similar to that derived in the 2XMM for 
HD 168112, another non-thermal radio emitter, at 
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its lowest luminosity (|Nazel l2009l) . CygOB2#9 
is thus only slightly harder and slightly more lu- 
minous than "normal" stars. However, it must 
be kept in mind that (1) with its long period, 
the Cyg 0B2 #9 binary is wide, and therefore the 
winds are quite diluted before colliding ; (2) the 
observations were taken rather far away from the 
periastron passage, at phases </)=0. 14-0.22 (Fig. 
[2|). Additional observations during such an event 
are needed before one can totally exclude a signif- 
icant contribution from wind-wind collision to the 
X-ray emission of Cyg 0B2 ^^9. 

Finally, it must be noted that CygOB2#9 
presents some variations between our observa- 
tions. While the overall spectral shape changes 
only slightly, the flux clearly decreased in the last 
two observations by about 10%, and there seems 
to be sonic shorter-term variations of lower am- 
plitude between the first four observations. The 
origin of such variations is unknown. The rea- 
son for the long-term 10% change cannot be con- 
strained before securing periastron observations, 
but it could simply be linked to a brightening as- 
sociated with strong wind-wind collisions around 
periastron (cf. the cases of WR140 and 77 Car). 

5. Conclusion 

Cyg 0B2 ^% is a rare case of non-thermal (ra- 
dio) emission associated with 0-stars. Such emis- 
sion is now thought to be associated with wind- 
wind collision in a binary, and we showed two years 
ago that Cyg 0B2 #9 was indeed a multiple sys- 
tem. However, understanding non-thermal emis- 
sion requires modelling the winds in detail, which 
in turn requests more than a "simple" binarity de- 
tection. 

The continuous monitoring of Cyg 0B2 #9 has 
led us to the derivation of a first orbital solution. 
The period is long, 2.4yrs, in agreement with the 
observed long-term radio modulation; the system 
eccentricity is large, 0.7-0.75, while the mass ra- 
tio approaches unity: Cyg 0B2 #9 is thus also 
one of the few known long-period 0-l-OB bina- 
ries and one of the few 0+OB systems presenting 
a high eccentricity. It should be noted that the 
RV curve of CygOB2#9 is peculiar, with only 
one unblended configuration (blueshifted primary 
- redshifted secondary) seen at periastron. Near 
apastron, the small RV difference between the two 



components is compatible with the large eccen- 
tricity and the orientation of the orbit, but the 
absence of RV crossing (or a limited one) requires 
in addition a large difference in the center-of-mass 
velocities for the two components. 

An additional monitoring was performed at 
high energies, with the hope of finding a signa- 
ture of a wind-wind collision (which should be 
the origin of the non-thermal radio emission). In 
the X-ray range, however, Cyg 0B2 #9 displays 
no large overluminosity nor any strong enhance- 
ment of its hard emission. There is thus, at least 
outside periastron, no clear, unquestionable signa- 
ture of X-ray emission from the wind-wind colli- 
sion. However, the flux varies, on both short and 
longer-timescales, with a 10% brightness decrease 
between 2004 and 2007. The cause of these varia- 
tions needs to be be investigated, notably by get- 
ting data closer to periastron. 

The next periastron passage of Cyg 0B2 #9 
should occur in June- July 2011, this time without 
any solar conjunction problem. It is the best time 
for finishing the study of this object by performing 
a multi-wavelength campaign (simultaneous radio. 
X-rays, and optical monitoring), which will finally 
open the possibility of modelling the rare high- 
energy phenomena occurring in this system. 
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Fig. 1.— Left: The He I A 5876 line in the Sophie 
data of HJD 4463.748 (thick sohd black line), 
4906.675 (thin solid red line, data taken half an 
orbit later than the previous ones), and 5142.240 
(dotted thick green line, data taken a quarter of an 
orbit later than previous ones) , as well as the Asi- 
ago data of the 2006 periastron passage (4022.252, 
dashed blue line) . Quoted phases are from the SB2 
solution (2nd column of Table [3]). Right: Deblend- 
ing using the method for these spectra. 
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Fig. 2. — Evolution of the Radial Velocities of 
He I A 5876 measured by the method with time 
(bottom panel) and phase (top pa nel, using the 
ephemeris from lvan Loo et al ]l2008l) . The primary 
(resp. second.) RVs are shown in filled (rcsp. 
open) circles ; vertical lines indicate 2006 Octo- 
ber 13 {HJD 4022.252), the approximate date 
of the periastron passage, and the dates 2.355yrs 
before/after. Dotted lines indicate the phases 
or dates at which the yLMM.- Newton data were 
taken. 
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Fig. 3. — RV curves superimposed on the best 
SB2 solutions (SB2 shown on the left panel, SB2' 
on the right panel). The symbols are as before 
and the errors on the velocities are of the order of 
10-20 kms-i. 



Cyg 0B2 #9; 29 October 2004 
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Fig. 4.— Left: EPIC spectra of Cyg 0B2 #9 and 
best-fit model in the first observation (lower lines: 
EPIC-MOSl in black + M0S2 in red ; upper line: 
pn in blue). Right: Variations of the fitted param- 
eters through the observations. 
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Table 1: Journal of observations. Heliocentric Ju- 
lian dates (mean values if N^^l) are given in the 
format HJD-2 450000, N is the number of spec- 
tra taken, AA is the wavelength range, R is the 
resolving power {X/ FWHMcaiib), S/N is the av- 
erage signal-to-noise ratio of the individual expo- 
sures around 5835A (or 4555A for the blue spec- 
trum). 



Instrument 


Date 


N 


AA (A) 


R 


S/N 


Espresso 


3288.806 


1 


3780-6950 


18000 


100 


Aurelie 


4560.563 


3 


5500-5900 


8800 


40 




4591.557 


3 


5500-5900 


8800 


80 




4620.472 


3 


5500-5900 


8800 


70 




4652.481 


3 


5500-5900 


8800 


70 




4680.444 


4 


5500-5900 


8800 


60 




4711.439 


1 


4450-4900 


7000 


30 




4715.388 


3 


5500-5900 


8800 


75 


WIRO 


4672.764 


4 


5300-6700 


4000 


200 




4746.743 


1 


5300-6700 


4000 


150 




4818.624 


1 


5300-6700 


4000 


220 




4842.070 


2 


5300-6700 


4000 


150 


Sophie 


4906.675 


3 


3900-6900 


35000 


65 




4936.968 


3 


3900-6900 


35000 


80 




4953.595 


1 


3900-6900 


35000 


75 




5006.497 


2 


3900-6900 


35000 


80 




5022.578 


2 


3900-6900 


35000 


80 




5058.433 


1 


3900-6900 


35000 


75 




5079.965 


2 


3900-6900 


35000 


70 




5123.370 


1 


3900-6900 


35000 


75 




5142.240 


1 


3900-6900 


35000 


95 
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Tabic 2: RVs of the He i line at 5875.62A, with he- 
liocentric Julian dates given in the format HJD- 
2450000. The RVs have been corrected for re- 
maining shifts using the DIB at 5780A (see text). 
Phases correspond to the SB2 solution (2nd col- 
umn of Table El). 



Date 


4>SB2 


i?V"i(kms-^) 


RV2(kms-'-) 


i?VGauss(kmS"l) 


3288.806 


0.12 


-60.3 


2.5 


-41.4 


3726.231 


0.63 


-52.6 


15.9 


-22.6 


3887.493 


0.82 


-47 


34 


-7. 


3990.455 


0.94 


-78.1 


70.2 




4022.252 


0.98 


-142.5 


96.1 




4036.237 


1.00 


-119.9 


130.9 




4051.389 


0.02 


-79.2 


112.7 




4066.298 


0.03 


-108 


19.6 




4244.476 


0.24 


-24.2 


-4.9 


-16.7 


4303.495 


0.31 


-27.9 


1 


-21.3 


4324.427 


0.34 


-36 


-0.4 


-19.6 


4348.765 


0.37 


-26.6 


-0.4 


-17.3 


4379.063 


0.40 


-25.9 


-10.1 


-17.1 


4463.748 


0.50 


-22 


1.8 


-10.1 


4472.247 


0.51 


-24.7 


-1 


-17.8 


4560.563 


0.61 


-29.1 


26.2 


-18.1 


4591.557 


0.65 


-18.8 


23.8 


8. 


4620.472 


0.68 


-23.9 


23.6 


14.2 


4652.481 


0.72 


-21.8 


19.8 


-1. 


4672.764 


0.75 


-26.3 


11.5 


-2.7 


4680.444 


0.76 


-22.4 


29.8 


13. 


4715.388 


0.80 


-11.4 


23 


1.3 


4746.743 


0.83 


-17.5 


24.2 


0.7 


4818.624 


0.92 


-30.8 


52.4 




4842.070 


0.95 


-48.4 


71.1 




4906.675 


0.02 


-97.1 


64.5 




4936.968 


0.06 


-68.9 


17.3 




4953.595 


0.08 


-54.7 


7.5 




5006.497 


0.14 


-40.3 


3.8 


-24.3 


5022.578 


0.16 


-36.5 


-10.9 


-25.9 


5058.433 


0.20 


-27.4 


-8.2 


-15.9 


5079.965 


0.22 


-32.7 


-6.8 


-17.2 


5123.370 


0.28 


-26.3 


-15.2 


-20. 


5142.240 


0.30 


-28.2 


-10.7 


-19.1 
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Table 3: Preliminary orbital solution for CygOB2#9. SBl' and SB2' refer to the use of identical velocities when the lines are blended (i 
text). 



Parameter 


SB2 


SBl for Prim. 


SBl for Sec. 


SBl for {RVi - RV2) 


SBl' for (RVi - RV2) 


SB2' 


P(d) 


852.9±4.3 


858.6±7.0 


848.4±3.0 


851.4±5.5 


851.3±5.6 


852.8±4.4 


u 


4036. 8±3. 6 


4019. 3±5. 7 


4047.6±2.3 


4036.5±4.4 


4031. 0±4. 9 


4030. 9±3. 9 


e 


0.744±0.030 


0.752±0.033 


0.799±0.033 


0.736±0.036 


0.704±0.034 


0.708±0.027 




-164.4±4.1 


167.1±5.8 


33.5±4.8 


-166.5±5.3 


-176.4±5.9 


-175.1±4.4 


M1/M2 


1.17±0.22 






1.27±0.35 


1.16±0.26 


1.10±0.17 


Ti (kms"""^) 


-40.6±3.2 


-38.5±2.0 




-38.9±4.0 


-24.4±5.3 


-28.1±3.4 


72(kms~^) 


16.6±3.5 




18.9±2.0 


17.8±4.0 


2.2±5.3 


1.2±3.5 


Xi(kms-i) 


53.0±7.0 


59.6±6.5 




50.3±11.0 


62.3±10.7 


63.2±6.8 


X2(kms-i) 


62.1±8.1 




73.2±7.7 


63.9±11.0 


72.0±10.7 


69.4±7.4 


ai sini(R0) 


598.0±84.4 






572.9±129.6 


744.3±132.7 


752.1±86.0 


02 sini(R0) 


699.3±97.7 






727.8±132.2 


860.2±134.3 


825.9±93.7 


Misin^i(MQ) 


21.7±7.2 






22.8±10.2 


41.0±15.4 


38.0±10.1 


M2sin^j(M0) 


18.6±6.1 






17.9±8.7 


35.5±14.0 


34.6±9.1 


rms(kms~^) 


14.5 


11.6 


8.1 


14.0 


12.5 


16.3 



Table 4: Journal of the 'K.MM-Newton observa- 
tions, with heliocentric Julian dates given in the 
format HJD-2 450 000 and phases from the SB2 
solution (2nd column of Table [3]). Xspec count 
rates (in cts~^) within the extraction region are 
given for the three instruments in the 0.4-10. keV 
energy band. 



Obs 


Date 




MOSl 


M0S2 


pn 


1 


3308.583 


0.14 


0.148±0.004 


0.162±0.004 


0.406±0.007 


2 


3318.558 


0.16 


0.138±0.004 


0.151±0.004 


0.397±0.008 


3 


3328.543 


0.17 


0.135±0.004 


0.141±0.004 


0.303±0.005 


4 


3338.505 


0.18 


0.131±0.004 


0.143±0.005 




5 


4220.355 


0.21 


0.108±0.004 


0.100±0.004 


0.301±0.007 


6 


4224.170 


0.22 


0.108±0.004 


0.109±0.003 





Table 5: Best- fit parameters models for 
CygOB2#9 and X-ray fluxes at Earth in 
the 0.5-10. keV band. 



Obs log Nwind 


kTi 


normi 


kT2 


norm2 xt (d.o.f.) 


fobs 


rcorr 


cm ^ 


keV 


lO^^cni-s 


keV 


10-3 cm-5 


ergcm"^ 


erg cm"^ s~^ 



1 


21.661:?^ 


0.671 


.12 
.04 


(3.4il°:«?) 


'?q+-27 


(1.891:1) 


0.86 (391) 


1.66 X 10" 


■12 


5.35 X 10- 


-12 


2 


21.561:11 


0.631 


.05 
.04 


(2.58l°-I°) 


2.2il:i^3 


(2.101:1) 


0.89 (423) 


1.60 X 10" 


12 


5.24 X 10- 


-12 


3 


21.69l:°6 


0.621 


.03 
.03 


(4.2ll°:f,) 


2-401:20 


(1.951:1) 


1.16 (414) 


1.74 X 10" 


12 


5.91 X 10" 


-12 


4 


21-80t:i3 


0.611 


.07 
.12 


(4.491?°^) 


2.69l:^? 


(1.801:^4) 


0.95 (122) 


1.69 X 10" 


12 


5.04 X 10" 


-12 


5 


21.85l:[J^ 


0.501 


.11 
.05 


(5.45li:fg) 


2.431:1 


(i.92i:i«) 


0.94 (211) 


1.52 X 10" 


12 


4.73 X 10" 


-12 


6 


2i.7ot:i!^ 


0.641 


.08 
.05 


(3.iolJ:°^) 


2 44+-34 


(i.8ol:^D 


0.79 (155) 


1.53 X 10" 


12 


4.67 X 10" 


-12 
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